Variational and Potential Formulation for Stochastic Partial Differential Equations 
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There is recent interest in finding a potential formulation for Stochastic Partial Differential Equa- 
tions (SPDEs). The rationale behind this idea lies in obtaining all the dynamical information of 
the system under study from one single expression. In this Letter we formally provide a general 
Lagrangian formalism for SPDEs using the Hojman et al. method. We show that it is possible to 
write the corresponding effective potential starting from an s-equivalent Lagrangean, and that this 
potential is able to reproduce all the dynamics of the system, once a special differential operator 
has been applied. This procedure can be used to study the complete time evolution and spatial 
inhomogeneities of the system under consideration, and is also suitable for the statistical mechanics 
description of the problem. 
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The inverse problem of variational calculus have 
been used by Hojman et al. |f( in the study of systems 
associated with first and second order deterministic 
differential equations. The propoused method permits 
to obtain all the dynamical infomation of the system 
allowing the quantization in terms of conserved quan- 
tities prescibed for the differential equation. On the 
other hand, a variational formalism has been devised 
by Gambar and Markus (see for example 0, 0, 3) 
as groundwork for proposing a field theory for non- 
equilibrium thermodynamical systems, giving valuable 
information about the entropy in terms of current den- 
sity and thermodynamic forces. Effective actions can 
be found by means of the Martin-Siggia-Rose formal- 
ism ||, a perturbative procedure that makes use of 
both physical and "conjugate" (auxiliary) fields. Also, 
Hochberg et al. have proposed in a series of papers 
0, 0, HI a "direct" approach, finding effective actions 
and potentials for SPDEs using a functional integral 
formalism similar in structure to those of quantum 
field theory. Last year Ao Q has reported a worth con- 
sulting novel approach for constructing potentials asso- 
ciated with first order SPDEs. All this efforts conduct 
to the idea that the variational formalism, in its direct 
or inverse form, could be the mathematical mechanism 
with the necessary tools for exploring the inherent dy- 
namics of physical, biological and chemical systems. 
In particular, structural development in cosmology 
and biology, pattern-formation, symmetry breaking, 
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population dynamics, chemical chaos and turbulence 
\WL Ull Ha . H3L II4L llq , are a few cases of a wide spec- 
trum of phenomena for which self-regulation, oscilla- 
tion, adaptiveness and multi-equilibrium behaviours 
can be adequately modeled by means of SPDEs. In a 
more recent paper |16| Hochberg et al. have reported 
that their effective potential does not provide infor- 
mation about the time evolution and spatial inhomo- 
geneities of the system under consideration, and the 
effective action is very difficult to calculate exactly. 

We will use in this Letter the Hojman et al. pro- 
ceedure, which enables to construct the Lagrangean 
for any regular mechanical system as a linear combi- 
nation of its own equations of motion. This particular 
construction is much wider than the traditional defini- 
tion L = T — V, which is only true when the "forces" 
involved are derivable from position-dependent po- 
tentials, therefore it may be used for general non- 
conservative systems. Its application to the study of 
SPDEs may provide additional understanding of the 
internal structure of these phenomena and also en- 
ables the use of a well known mathematical machinery 
to find conserved quantities, equilibria and stability 
cases, and other dynamical properties. For an exam- 
ple of the application of this method to self-regulated 
systems, refer to |l7j . 

The goal of this Letter is to establish a varia- 
tional and effective potential formalisms for the study 
of SPDEs with arbitrary noise function. This ap- 
proach have the particular advantages that it provides 
Lagrange-Hamilton functionals in a very direct way 
starting from the equations of motion, and that the 
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effective potential deduced is able to take account of 
the stability and equilibria conditions, temporal evolu- 
tion, anisotropies and inhomogeneities of the system. 

We will consider in this work SPDEs that can be 
written as 



Uip 1 - $ l o j ') -C =0 



(1) 



for i.j = 1, ...,m where m is the number of degrees 
of freedom, i\f = ip l (qi , t) denotes the components 
of a vector field whose arguments are, in the general 
case, spatial coordinates q J and time t; □ is an ar- 
bitrary linear space or time (or both) differential op- 
erator that does not depend on the field ip l \ <fr l (^ J ) 
is any, ussually non-linear, deterministic forcing term 
and = (q J , t) is a random function of its arguments 
describing the stochastic force (noise) in the system. In 
Table 1 some particular cases of □ operators and $ l 
functions are presented (a p are the coefficients of the 
polynomial of order p while k and u> are real scalars; 
for details see and references therein). Hereafter, 
Einstein summation convention and Euclidian metric 
tensor are assumed. The noise functions to be consid- 
ered here are completely arbitrary. 

It is important to bear in mind that even when 
it is true that different sorts of noise and differ- 
ent behaviours of the underlying deterministic par- 
tial differential equations (hyperbolic, parabolic, lin- 
ear or nonlinear) require different techniques to find 
the corresponding solutions, those may be unneces- 
sary, as we will presently show, for obtaining a varia- 
tional/effective potential formulation of the problem. 



Operator: 


□ 


Function: 




D'Alembertian 


9?-V 2 


Polynomial 




Diffusion 


d t - kV 2 


Burgers (noisy) 


<§ (v^) 2 


Temporal 


dt 


Purely Dissipative 


Sip' 



Table 1. Some Operators and Dissipation Functions. 

The equation of motion for a mechanical system 
arises from a set of m differential equations. For con- 
ceptual reasons, it is suitable to rewrite as equa- 
tions of motion in the variational sense 



i> l -F l (y,9 fc V> j ,V>V) =0 



(2) 



where F % behaves as "forces" (both deterministic and 
stochastic) divided by unitary mass and may include 
spatial derivatives of the field; the dot stands for total 
temporal derivative. 

In the Inverse Problem of the Variational Calcu- 
lus the Lagrangean , ■0 J , t) is constructed such 



that relations @ can be effectively deduced via Euler- 
Lagrange equations. The existence of such a La- 
grangean is studied using the Helmholtz conditions 0] 
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Nevertheless, these conditions do not give any war- 
ranty about uniqueness. Two Lagrangeans are said to 
be solution-equivalent (or s-equivalent) if they differ 
only by a global multiplicative constant, 77, and a total 

time derivative of some gauge 



i]L = L 



~dt 



(4) 



The different systems of equations they provide, 
however, have exactly the same equations of motion. 

The Hojman et al. method enables us to write L as a 
linear combination of the known equations of motion; 
then for i, j = 1, to, 



where 
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(5) 
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(7) 



In equation © the quantities under partial derivative 
(fl 2m ) are constants of motion of the system, while the 
corresponding coefficients (2?2m— 1) are arbitrary func- 
tions whose arguments are constants of motion. There 
are plenty of ways to write these £*2m functions. For 
instance, one possible form for the T>i m -\ functions, 
given the D? ™, conserved quantities, is presented in 
reference |19). 

When the conserved quantities are unknown the 
problem is reduced to find ^ such that the following 
system 



dip 3 dip 



di 



dF 3 

dipi 



is satisfied, and 
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^0 (9) 
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where the on-shell derivative ^ is defined such that it 
behaves as an usual total time derivative taking into 
consideration that, on the shell, equation Q is always 
true. 

It is important to remark that this method is useful 
for both second-order and first-order differential equa- 
tions [D. For further details the reader is exhorted to 
review 0, 0] and the references therein. 

By virtue of equation JSJ, the general Lagrangean 
for SPDEs of the form JIJ can be written as 



L = m {Utf - $ l - C 



(10) 



where the [ii parameters must be determined for each 
case in study. 

It can be shown that the 2m constants of motion of 
equation IjlUI) can be formally cast in 



D {2i-1) = ^i_ / ^i dt 



F l dt 



(11) 



(12) 



It is clear that equations (|llf> and l|12l) must satisfy 
conditions 10)- © of the Hojman et al. method. 

Once the Lagrangian I|1U|I is completely determined, 
the corresponding Hamiltonian can be found trivially 
by the usual Legendre transformation, and, as will be 
shown in the following lines, also the effective potential 
can be written straightforwardly. 

Let us define now the following differential operator 



* 1 d 
Vi = 

Ipi dt 



(13) 



Then, let require that the effective potential, V e ff, 
be such that 



i v eff 



-Ft 



(14) 



Note that in this case the on-shell derivative takes the 
form 



d d d , d 
— = h ib 1- F — - 

dt dt v di>* 



(15) 



and thus, in the conservative case, V e ff = V e ff(ip) = 
V, equation l|14|) provides (from now on the symbol 
= indicates that the relation is valid only in special 
cases) 



dV 

ViV e ff ■ — F 



dtp 1 



(16) 



which is equivalent to the standard relation T/iV — 
—Fi when the vector field ip l coincides with the spatial 
coordinate q l . 

Now, wc arc interested in finding a gauge A — X such 
that equations H1(J|) and (jSJ provide an s-equivalent 
Lagrangean, L, that can be written as the difference 
between some function (kinetic energy) T e ff and the 
effective potential V e ff- In consequence (for the sake 
of simplicity, r) = 1) 

Veff = T eff - w (D^ - & - f) + § (17) 

For a general non-conservative system, T e f / can be 
written as [20| 



T. 



f f = aA[f)+pB + 7 g(^,#) (18) 



and thus, it is necessary to determine the scalars a, j3, 7 

and the functions A (^) , B (ijj 1 , ^ , G (^tpi^A to 

completely define Veff- 

It should be useful to write the kinetic energy of the 
system in the traditional quadratic form 



1 



l eff 



(19) 



In order to do so, it is necessary to find the correspond- 
ing gauge first. Equation l|14|) provides the necessary 
constraint; thus, by taking the nabla-star derivative 
at both sides of equation (|17J) we obtain, for on-shell 
trajectories, 



V; 



dX 
~dt 



H :l (Uib ] - <1> ' - ^ 



-2K- 



(20) 



Solving this equation for the time derivative of the 
gauge A, we have 



= -2 / ftFidt + Mi OV ~ ~ D + V Q (21) 



dX 
~dl 



where Vo is an arbitrary constant of integration. Con- 
sequently, following equation l|17fl 



V e ff =V + ^ipi 



(22) 



which is the general effective potential for equation (JTJ 
given the choice (|19|l for the kinetic energy. 

In summary, we have presented in this Letter a novel 
way to provide both variational and effective potential 
formulations for general SPDEs with arbitrary noise 
function. There are several useful applications of this 
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result. Once the Hamiltonian is obtained, for exam- 
ple, quantization of systems described by equation |T|1 
follows straightforwardly. Also, the Hamilton- Jacobi 
approach may help to solve the equations of motion of 
a system via a convenient, if possible, variable separa- 
tion. 

The effective potential 11170 is constructed such that 
it contains all the dynamical information of the sys- 
tem. In general, it may have explicit dependence on 
time, on the field itself and on the derivatives of the 
field; certainly, the nabla-star operator identifies the 
contribution of each functional dependence by means 
of specific terms, as can be seen in equation <|15[) . 

As reported 0, 0] , in other approaches the effec- 
tive potential is useful only for stationary or static 
regimes; the approach presented here can be used to 
study those cases, the complete temporal evolution of 
and also the equilibria and stability states. The 
present potential formulation is based on a classical 
mechanics approach, there is no need of auxiliary or 



ghost fields. Nor Fokker-Planck equations neither spe- 
cial assumptions about the noise function are invoked 
in the construction. 

The important relation between the effective po- 
tential (|17l) and the statistical mechanics deserves a 
detailed discussion, and it will be treated elsewhere. 
However, it can be shown that this potential appears in 
the steady state solution of the corresponding Fokker- 
Planck equation (or even the Klein-Kramers equation) 
[jj , a Boltzmann-Gibbs distribution 

p (V' 4 ) = |cxp(-T/ e// (^)) (23) 

with Z as the partition function. 

A detailed application of this procedure to several 
particular cases of SPDEs, such as Langevin, inhomo- 
geneous Klein-Gordon, and diffusion equations from a 
mechanical, statistical and quantum point of view is 
in course of preparation. 
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